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Recently experimental techniques, such as magnetic force 
microscopy (MFM), have enabled the magnetic state of in- 
dividual sub-micron particles to be resolved. Motivated by 
these experimental developments, we use Monte Carlo simu- 
lations of two-dimensional kinetic Ising ferromagnets to study 
the magnetic relaxation in a negative applied field of a grain 
with an initial magnetization mo = +1. The magneto- 
static dipole-dipole interactions are treated to lowest order by 
adding to the Hamiltonian a term proportional to the square 
of the magnetization. We use droplet theory to predict the 
functional forms for some quantities which can be observed by 
MFM. One such quantity is the probability that the magneti- 
zation is positive, which is a function of time, field, grain size, 
and grain dimensionality. The relaxation is characterized by 
the number of droplets larger than a field-dependent critical 
size which form during the switching process. Our simula- 
tions of the kinetic Ising model are in excellent agreement 
with droplet-theoretical predictions. The qualitative agree- 
ment between experiments and our simulations of switching 
in individual single-domain ferromagnets indicates that the 
switching mechanism in such particles may involve local nucle- 
ation and subsequent growth of droplets of the stable phase. 

I. INTRODUCTION 

The processes by which magnetization reversal occurs 
in the nanoscale ferromagnets that will make up the next- 
generation recording media are the subject of active re- 
search. One quantity for which theory and experiment 
often disagree is the lifetime r, which is the time required 
for a particle with initial magnetization mo = +l to reach 
777. = when a magnetic field in the — z direction is ap- 
plied. Micromagnetics, Q a theoretic technique in which 
differential equations are numerically solved on a coarse- 
grained lattice, predicts that the lifetime is given by the 
Arrhcnius equation 

r oc exp(/3AF) (1) 

with AFocL d . Here is the temperature in units of 
energy, AF is the free-energy barrier between the sta- 
ble and metastable phases, and L is the linear system 
size. This same prediction is made by the standard Neel- 
Brown theory of single-domain ferromagnets. j^,|3| The 
evident failure of Eq. (|l|) with AF cx L d for somewhat 
larger grains is ascribed to the existence of more than one 
domain in larger particles, as is a corresponding peak in 



plots vs. L of the switching field H sw , which is the field 
required to yield a given lifetime. 

Recently techniques such as MFM have been used 
to resolve the magnetic properties of isolated, well- 
characterized single-domain particles (see, e.g., Ref. 
This is an important advance, since previous experiments 
on ferromagnetic powders left uncertainties due to the 
range of grain sizes and orientations and the local mag- 
netic environments. Observations of individual particles 
by MFM have made it clear that even for many single- 
domain particles, Neel-Brown theory is inadequate. 

We have applied the statistical-mechanical droplet the- 
ory of metastable decay to nanoscale ferromagnets with 
large uniaxial anisotropy, and compared Monte Carlo 
simulations of square-lattice Ising systems with droplet- 
theory predictions. § (For a review of droplet theory, 
see Ref. ||.) The agreement between theory and simula- 
tion is quite good, and despite the crudeness of the Ising 
model as a model for real magnets, it shows good qualita- 
tive agreement with the MFM experiments. We find rich 
L-dependent behavior in the standard Ising model, even 
though its equilibrium structure is a single domain for 
all L. This suggests that for some strongly anisotropic 
magnetic materials, magnetization reversal may occur 
through the nucleation and growth of non-equilibrium 
droplets. Details of our work are given in Refs. § and 
0, including formulae for general dimensionality. For 
simplicity, we only discuss the two-dimensional case here. 

II. APPLIED DROPLET THEORY 

To be concrete, consider a two-dimensional kinetic 
Ising ferromagnet (s, = ±l) with Hamiltonian 

H = -jY,^- H J2 s * + L ~ 2D (l2 s i) ( 2 ) 

n.n. i ^ i ' 

and Metropolis single-spin-flip dynamics on a square lat- 
tice with periodic boundary conditions. The last term 
represents a mean-field approximation for the dipole- 
dipole interaction energy and is taken to be zero except 
as noted. The critical radius of a "droplet" of Sj = — 1 
spins surrounded by Si = +1 spins occurs when the free 
energy of the droplet (2irRa — irR 2 2\H |) is maximum: 
R c w a/2\H\, where a is the surface tension per unit 
length. Droplets smaller than this will very probably 
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shrink and vanish; larger droplets will very probably grow 
and reverse the magnetization of the system. In a suffi- 
ciently large system, the probability per unit time that 
a critical droplet forms, centered at a given site, is given 
by droplet theory as M 



Ioc \H\ 3 exp(-/3ir<T 2 /2\H\) 



(3) 



The details of the magnetization reversal depends on the 
number of critical droplets the system forms. 
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FIG. 1. The relation between the switching field H aw 
and the system width L for two different fixed lifetimes 
(solid curves), calculated by kinetic Ising model simulations 
at k B T = 0.8fc B T c ~ 1.815J. The dotted curve is near the 
crossover between the CE and SD regions; || the dash-dotted 
curve is near the crossover between the SD and MD regions. 



For weak fields or small systems (L < B, c ), no critical 
droplet can form. This is called the Coexistence (CE) 
Region, and since two interfaces (remember: periodic 
boundary conditions) must form to reverse the magne- 
tization, 



exp {(3 [2a L - O (HL 2 )]} 



(4) 



For slightly larger L, the first supercritical droplet will 
grow to the size of the system before another one can 
form. The lifetime in this Single Droplet (SD) Region is 



(5) 



In both the CE and SD regions, switching is a Poisson 
process, so the standard deviation of the lifetime is com- 
parable to r. Both Eq. (ji|) and Eq. (||) are actually special 
cases of the Arrhenius equation, Eq. (|l|), but in neither 
case is AF proportional to L d . Note that if r is held 
constant and the system size is increased, Eq. (||) implies 
that the magnetic field must increase, whereas Eq. (||) im- 
plies that the magnetic field must decrease. This shows 
that the peak in H sw occurs near the crossover between 
the CE and SD regions (see Fig. Ill) . 




FIG. 2. The probability that m > for a kinetic Ising 
system in the SD region. T = 0.8T C , t = 914 Monte Carlo 
steps per spin (MCSS) and L — 10. The solid curve is a fit of 
exp(— t/r) to the MC data, where r is given by Eq. (§). The 
inset figure shows the fitted curve over a wider range in H. 
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FIG. 3. The probability that m > for a kinetic Ising 
system in the MD region. T = 0.8T C , r = 40.7 MCSS and 
L = 30, 100, and 300. The solid curves are fits of droplet 
theory predictions to the MC data. The dashed curve is the 
fit of the droplet-theory prediction for L — 100 extrapolated 
to L = 300. 

The probability that the magnetization is greater than 
zero P(m>0) is shown as a function of field in Fig. |^ for 
a system in the SD region. This probability is what is 
most easily observed in MFM experiments, and it decays 
exponentially with time in both the CE and SD regions. 
In the SD region, the system is very unlikely to return to 
the the metastable state from the stable state, so 

P(m>0) = exp (-t/r) . (6) 

In the CE region such backwards switching takes place 
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on a timescale comparable with the initial decay, so the 
situation is more complicated. 

For sufficiently large L or H, several supercritical 
droplets may form before any one of them has grown to 
the size of the system. This is the Multi-Droplet Region 
(MD). Such systems were first studied by Kolmogorov, 
H Johnson and Mehl, || and Avrami, [jlO| and have a 
lifetime 

1/3 



[/7rw 2 /31n2] 



(7) 



where the radial growth velocity v is assumed to be pro- 
portional to \H\. Although r is independent of L, the 
variance of the lifetime is proportional to (v/L) 2 . Mea- 
suring P(m > 0) as a function of H or t thus provides 
a means of estimating the proportionality constant be- 
tween v and H (see Fig. ||). Details are given in Ref. ||. 

The addition of the dipole-dipole interaction energy in 
Eq. (Q) makes the forms for the lifetime somewhat more 
complicated than we have presented here. In the CE and 
SD regions, a form of the Arrhenius equation [Eq. ([I])] 
still applies, and switching is still a Poisson process, so 
Eq. (^|) still applies for the SD region. In the MD region, 
the system evolves in a time-dependent effective field, 
H e g = H — 2Dm(t). The time-dependent magnetization 
can be found to O(D) fairly easily, and we find numer- 
ically that the 0(D 2 ) correction is relatively small. We 
can then solve analytically for r to 0(D 2 ) and find good 
agreement with simulation results (Fig. |]). A detailed 
treatment of the D > case is given in Ref. M . 
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explained by droplet theory and show good qualitative 
agreement with experiments, despite the comparative 
simplicity of the Ising model. This simplicity, in turn, 
allows us to develop an understanding of the underly- 
ing statistical mechanics. Particular features to make 
the model more realistic, such as appropriate bound- 
ary conditions, quenched randomness, and less rigorous 
anisotropy will be added in later studies. 
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FIG. 4. The lifetime r for a kinetic Ising system in the 
MD region, normalized by the lifetime in a similar system 
with D = 0. T = 0.8T C . The solid curves are droplet-theory 
predictions. 



III. CONCLUSION 



We have used Monte Carlo methods to simulate mag- 
netization switching in two-dimensional kinetic Ising fer- 
romagnets. The results of the simulations can be well 
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